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1. Introduction 

Color confinement is still a puzzling problem not withstanding the large mess of numerical 
investigations aimed to understand the nature of the QCD vacuum. Indeed, the mechanism 
that leads to color confinement remains an open question despite intense lattice studies for 
nearly three decades. 

According to a model conjectured long time ago by G. 't Hooft [1] and S. Mandel- 
stam [2] the confining vacuum behaves as a coherent state of color magnetic monopoles, or, 
equivalently, the vacuum resembles a magnetic (dual) superconductor. Up to now there is 
numerical evidence [3-12] in favor of chromoelectric flux tubes in pure lattice gauge vacuum. 
As well there have been extensive numerical studies [13-24] of monopole condensation. 

An alternative model for color confinement is based on the special role of center vor- 
tices and Z]\f symmetry (see [25] and references therein), even if the connection of center 
symmetry to confinement has been recently questioned in [26,27]. 

One may conclude that there is no totally convincing explanation of the confinement 
phenomenon (for recent reviews on confinement see [25, 28, 29]) and that a full under- 
standing of the QCD vacuum dynamics is still lacking. Indeed, as recently observed [30] in 
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connection with dual superconductivity picture, even if magnetic monopoles do condense in 
the confinement mode, the actual mechanism of confinement could depend on additional 
dynamical forces. Therefore we feel that it is important to explore any new paths that 
possibly may suggest new hints for understanding the QCD vacuum. 

In a previous paper [31] we reported numerical results showing that in four dimen- 
sional U(l) lattice gauge theory the confining vacuum behaves as a coherent condensate 
of Dirac magnetic monopoles, according to analytical results in the literature [32]. In the 
same paper we gave account of numerical results indicating condensation of abelian mag- 
netic monopoles and abelian vortices in the confined phase of finite temperature SU(2) and 
SU(3) lattice gauge theories in (3-1-1) dimensions. Therefore one might conclude that in 
SU(2) and SU(3) gaTigc theories the confining vacuum behaves as a coherent abelian mag- 
netic condensate. Wc also found [33] that a weak constant abelian chromomagnetic field 
at zero temperature is completely screened in the continuum limit, while at finite temper- 
ature [34, 35] our numerical results indicate that the applied field is restored by increasing 
the temperature. These results strongly suggested that the confinement dynamics could 
be intimately related to abelian chromomagnetic gauge configurations. Similar arguments 
have been reported in ref. [36]. Moreover, in Refs. [34,35] the SU(3) vacuum was probed by 
means of an external constant abelian chromomagnetic field with increasing field strength. 
Remarkably, wc found that by increasing the strength of the applied external field the de- 
confinement temperature decreases towards zero. This means that strong enough abelian 
chromomagnetic fields destroy the confinement of color. In analogy to what happens in 
familiar superconductors when the strength of an external magnetic field is increased (see 
for instance ref. [37]), this effect can be named "reversible color Meissner effect" . Altough 
the existence of a critical chromomagnetic field is not easily understandable within the co- 
herent magnetic monopole condensate picture of the confining vacuum, it could be directly 
explained if the vacuum behaves as an ordinary relativistic color superconductor, or differ- 
ently stated, if the confining vacuum resembles as a coherent condensate of tachionic color 
charged scalar fields. Thus we have to reconcile two apparently different aspects. Prom one 
hand, the confining vacuum docs display condensation of both abelian magnetic monopoles 
and vortices, on the other hand the relation between the deconfinement temperature and 
the applied abelian chromomagnetic field suggested that the vacuum behaves as a conden- 
sate of an effective color charged scalar field whose mass is proportional to the inverse of 
the magnetic length [34,35]. The reversible color Meissner effect could be in agreement 
with R. P. Feynman who, in a seminal paper [38], argued that in three dimensional SU(2) 
gauge theory long range correlation between gluonic degrees of freedom destroys confine- 
ment. We would like to point out, to avoid misunderstanding, that our reversible color 
Meissner effect is not related to color superconductivity in cold dense quark matter (for a 
recent review see ref. [39] and references therein). 

The aim of the present paper is to investigate if our reversible color Meissner effect is 
a generic feature of non abelian gauge theories. To this end, we shall compare SU(3) and 
SU(2) gauge theories in an external abelian chromomagnetic field both in (3+1) and (2+1) 
dimensions. Wc shall, also, consider three and four dimensional U(l) gauge theories in a 
magnetic background field. 
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The plan of the paper is as follows. In sect. § we briefly recall for reader convenience our 
proposal of lattice effective action and define the abelian chromomagnetic field on the lat- 
tice. In sect. ^ we present our results on vacuum dynamics in an external chromomagnetic 
background field in (3+l)-dimensions for SU(3) and SU(2) at finite temperature, and for 
U(l) at zero temperature. Sect. ^ is devoted to corresponding results in (2+l)-dimensions. 
Finally in sect. ^ we summarize and conclude. In Appendix |A| we present results for SU(3) 
in an external chromomagnetic background field directed along the direction 8 in color 
space. 

2. The gauge invariant lattice effective action 

In our previous studies, in order to investigate vacuum structure of lattice gauge theories 
both at zero and finite temperature, we introduced a lattice effective action for gauge 
systems in external static background fields. In this section, for reader convenience, we 
shall briefiy summarize our proposal of lattice effective action which is gauge invariant 
against static gauge transformations of the background field. 

2.1 The lattice eflfective action: T = 

In Refs. [33, 40] we introduced a lattice gauge invariant effective action r[A'^^*] for an 
external background field A^^^: 

where Lt is the lattice size in time direction and A^^^{x) is the continuum gauge potential 
of the external static background field. ^[A*^^*] is the lattice partition functional 

= [ VU e"^^ , (2.2) 

with Sw the standard pure gauge Wilson action. 

The functional integration is performed over the lattice links, but constraining the 
spatial links belonging to a given time slice (say Xf = 0) to be 

Uk{x,xt = 0) = Ut\x) , {k = 1,2,3) , (2.3) 

C/|^*(x) being the lattice version of the external continuum gauge potential A^^^{x) = 
A^^{x)Xa/2. Note that the temporal links are not constrained. 

In the case of a static background field which does not vanish at infinity we must also 
impose that, for each time slice xt ^ 0, spatial links exiting from sites belonging to the 
spatial boundaries are fixed according to eq. (p.3|). In the continuum this last condition 
amounts to the requirement that fluctuations over the background field vanish at infinity. 

The partition function deflned in eq. ( |2.2D is also known as lattice Schrodinger func- 
tional [41,42] and in the continuum corresponds to the Feynman kernel [43]. Note that, at 
variance with the usual formulation of the lattice Schrodinger functional [41,42] where a 
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lattice cylindrical geometry is adopted, our lattice has an hypertoroidal geometry so that 
Sw in eq. (^) is allowed to be the standard Wilson action. 

The lattice effective action r[^^^*] corresponds to the vacuum energy, £'o[j4'^^*], in 
presence of the background field with respect to the vacuum energy, -Eo[0], with A^^^ = 

p^^ext] — ^ So[i*''^*] -^o[0]. (2.4) 

The relation above is true by letting the temporal lattice size Lt ^ oo; on finite lattices 
this amounts to have Lf sufficiently large to single out the ground state contribution to the 
energy. 

Since the lattice effective action eq. (|2.lD is given in terms of the lattice Schrodinger 
functional, which is invariant for time-independent gauge transformation of the background 
field [41,42], it is gauge invariant too. 



2.2 The thermal partition functional 

If we now consider the gauge theory at finite temperature T = l/{aLt) in presence of an 
external background field, the relevant quantity turns out to be the free energy functional 
defined as 

^[1-1.4,4211^1. (2,5) 



1 ZtM 



ZtIA^^^] is the thermal partition functional [44] in presence of the background field A^^^, 
and is defined as 



Zt 



A 



ext 



PC/e""^^. (2.6) 



Ineq. B), as in eq. ( |2.2D , the spatial links belonging to the time slice xt = are constrained 



to the value of the external background field, the temporal links are not constrained. 
On a lattice with finite spatial extension we also usually impose that the links at the 
spatial boundaries are fixed according to boundary conditions eq. (p.3|), apart from the 
case in which the external background field vanishes at spatial infinity (as happens for the 
monopole field), where the choice of periodic boundary conditions in the spatial direction 
is equivalent to eq. (|2.3| ) in the thermodynamical limit. If the physical temperature is 
sent to zero, the thermal functional eq. reduces to the zero-temperature Schrodinger 
functional eq. ( |2.2| ). The free energy functional eq. (|2.5| ) corresponds to the free energy, 
^j^extj^ in presence of the external background field evaluated with respect to the free 
energy, F[0], with ^4*^^* = 0. When the physical temperature is sent to zero the free energy 
functional reduces to the vacuum energy functional eq. (|2.1|). 

2.3 Abelian (chromo) magnetic background field 

We are interested in vacuum dynamics of U(l), SU(2), and SU(3) lattice gauge theories 
under the infiuence of an abelian chromomagnetic background field. 

In our previous studies we found that in SU(2) and SU(3) at zero temperature a 
(not too strong) constant abelian chromomagnetic field at zero temperature is completely 
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screened in the continuum limit [33]. We also found that in SU(3) the deconfinement tem- 
perature depends on the strength of an applied external constant abelian chromomagnetic 
field [35]. This is at variance of abelian magnetic monopoles where the abelian monopole 
background fields do not modify the deconfinement temperature [45]. We would like to 
corroborate our findings with further investigations, in particular we would like to ascertain 
if the dependence of the deconfinement temperature on the strength of an applied external 
constant abelian chromomagnetic field is a peculiar feature of non abelian gauge theories. 

Let us now define a static constant abelian chromomagnetic field on the lattice. We 
first consider the SU(3) case. In the continuum the gauge potential giving rise to a static 
constant abelian chromomagnetic field directed along spatial direction 3 and direction a in 
the color space is given by 



(2.7) 



In SU(3) lattice gauge theory the constrained lattice links (see eq. (|2.3|)) corresponding to 
the continuum gauge potential eq. ( p.Tj ) are (choosing a = 3, i.e. abelian chromomagnetic 
field along direction 3 in color space) 



C/f^*(x) 
Ul^\x) 



X) = 1. 

gHxi\ 







exp(i 2 

exp(— i 2 




(2. 



We will refer to this case as T3 abelian chromomagnetic field. If we choose instead abelian 
chromomagnetic field along direction 8 in color space the constrained lattice links are given 

by 



Ul''\x) = Ul''\x) = 1 
9. 

2v/3 



uT\x) 



exp(i2^' 



exp(i 





gHx\ - 









2v^ ' 

exp(-z2ff: 



(2.9) 



We will refer to this case as Tg abelian chromomagnetic field. Since our lattice has the 
topology of a torus, the magnetic field turns out to be quantized 



ogH 27r 
a — = -^next 

Z hi 



next integer. 



(2.10) 



In the case of SU(2) lattice gauge theories the constrained spatial links are 



ur{x) = ur{x) = i, 

C/|-(^) = cos(^) + za=^sin(^), 



(2.11) 



cT'^ being the Pauli matrix. 

Finally in the U(l) case the constrained spatial links corresponding to a constant magnetic 
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background field (along spatial direction 3) are 

C/|^t(f) = cos{gHxi) + ism{gHxi) . 

Since the free energy functional is invariant for time independent gauge transforma- 

tions of the background field A^^^, it follows that for a constant background field, is 
proportional to the spatial volume V = Lf, and the relevant quantity is the density /[^'^^*] 
of free energy 

/[l^^t] = -^^[1^^*] . (2.13) 

We evaluate by numerical simulations the derivative with respect to the coupling (3 of the 
free energy density /[^°^*] at fixed external field strength gH 

^,[^ext]^ /l ^ lj^^Trt/^,(x)\ -/l iReTr[/,,(x)\ , (2.14) 

\ X,fJ.<U / Q \ X,fJ.<U I ^ext 

where the subscripts on the averages indicate the value of the external field and = L^x Lt 
is the lattice volume. The generic plaquette U^y{x) = Ufj,{x)Uy{x + fl)uji{x + i')ul{x) 
contributes to the sum in eq. (|1D if the hnk U^{x) is a " dynamical" one, i.e. it is not 



constrained in the functional integration eq. (|2.6| ). Observing that /[^d*^^*] = at /3 = 0, 
we may eventually obtain /[^d*^^*] from /'[A*^^*] by numerical integration: 

/[yl-^^t] = / . (2.15) 



JO 

3. (3+1) dimensions 

In this section we report results obtained in studying the finite temperature phase transition 
of lattice gauge theories SU(3) and SU(2) in (3-|-l)-dimensions, in presence of a constant 
abelian chromomagnetic background field. We shall also report results for confinement- 
Coulomb phase transition in U(l) lattice gauge theory at zero temperature in a constant 
magnetic background field. A preliminary account of our results has ben presented in 
ref. [46]. 

3.1 SU(3) 

We simulate SU(3) pure gauge theory in a constant abelian background field defined in 



Eqs. (2.7) and (2.8). As is well known, the pure SU(3) gauge system undergoes a decon- 
finement phase transition at a given critical temperature. Our aim is to study the possible 
dependence of the critical temperature from the strength of the applied field. The critical 
coupling (3c can be evaluated by measuring /'[A''^*], the derivative of the free energy den- 
sity with respect to /3, as a function of j3. Indeed we found that /'f^*^^*] (see eq. ( |2.14|) ) 
displays a peak in the critical region (see fig. ||) where it can be parameterized as 

/'(/3,Lt) ai{Lt) 



'ext 



a2{Lt)[(3-f}*{Lt)Y + l 



(3.1) 
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Figure 1: SU(3) in (3+1) dimensions. The derivative of the free energy density with respect to the 



gauge couphng /3, eq. ( |2.14| ), versus (3 at fixed external field strength (rjcxt — 1) for spatial lattice 
size Ls = 64 and temporal lattice sizes Lt — 4, 5, 6, 7, 8. Solid lines are the fits eq. (|3.l[). 



In eq. (|3.lD we normalize /' to Sg^t > the derivative of the classical energy due to the external 
applied field 



£ 



,2tt 



ext ^ 3 ~ — )] = 3 [1 - C0S( — next)] ■ 



(3.2) 



Remarkably, we have checked that the evaluation of the critical coupling f3*{Lt) by means of 
y/j^extj -g consistent with the usual determination obtained through the temporal Polyakov 
loop susceptibility: 



X{\P\)=<\P?>-<\P\ 
^ = ^E^Tr n C/4(x4,x). 



(3.3) 



2:4=1 



The Polyakov loop susceptibility near the peak has been obtained by means of the density 
spectral method [47,48]. The statistical errors for the points belonging to the extrapolated 
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Figure 2: SU(3) in (3+1) dimensions. The susceptibility of the absolute value of the Polyakov loop 
evaluated for three values of /3 on a 64'^ x 6 lattice and external field strength next = 1- The solid 
line has been obtained by reweighting. The error bars have been estimated using the bootstrap 
method. 



curve near the peak, as well as the position of the peak and its statistical error, were 
evaluated by means of a bootstrap analysis [49]. For instance, on a 64^ x 6 lattice and 
'^ext = 1 we get Pc = 5.6272(69) from eq. (|3.l|) and Pc = 5.6266(12) when evaluating 
the peak of the Polyakov loop susceptibility by means of the density spectral method (see 

%•!)• 

Once (3*{Lt) has been determined, the deconfinement temperature can be preliminarily 
estimated in units of Ai^tt- Indeed 



Tc 1 1 



Alatt Ltfsu{3){/3*{Lt))' 

where 



(3.4) 
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N being the color number, bo = (UN) / {AStt"^) , and h = {34:N^)/{3{l6TT^ f). In order 
to obtain the continuum hmit of the critical temperature we have to extrapolate Tc/ Aiatt ; 



given by eq. (3.4), to the continuum. This can be done, following ref. [50], by means of a 
linear extrapolation of T^/Aiatt as a function of aT^ for aT^ 0. We varied the strength 
of the applied external abelian chromomagnetic background field to study quantitatively 
the dependence of Tc on gH. From fig. ^, where we display Tc/Aiatt versus gH in lattice 
units, we may conclude that the critical temperature decreases by increasing the strength of 
the external abelian chromomagnetic field and eventually goes to zero for a strong enough 
external field. 

To get more insight into this result we can try to parameterize the behavior of the 
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Figure 4: SU(3) in (3+1) dimensions. Tlic critical temperature Tc estimated on a 64'^ x 8 lattice 



in units of the string tension, eq. (3.1C), versus the square root of the field strength \/gH in units 
of the string tension, eq. ( p.ll ). Solid line is the linear fit eq. (3.12). In correspondence of zero 
vertical axis: open circle is Tc/ y/a at zero external field, Eq. (3.13); full circle is the determination 
of Tc/ \fa obtained in the literature [51] The green full point is the critical field in units of the string 



tension, eq. (3.14) 



critical temperature versus the applied field strength. As a matter of fact, if the magnetic 
length, defined as an ~ 1/ ^fgH^ is the only relevant scale of the problem for dimensional 
reasons one expects that 

Tl ~ gR. (3.6) 
Indeed, as fig. ^ displays, we get a good fit to our data using the following parameterization 



Tc(gg) T,(0) 
— 7 = —. V ay a'^gtl . 



(3.7) 

^>-latt ^»-latt 
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We get 

7;(0)/Aiatt = 35.5 ± 5.2 

(3.8) 

a = -42.4 ±7.4. 

It is worthwhile to note that our estimation for Tc(0)/Aiatt is compatible with rc(0)/Aiatt = 
29.67 lb 5.47 obtained in ref. [50] with completely different methods. 

The preliminary analysis of our lattice data drives us to conclude that, remarkably, a 
critical field, gHc ~ 0.68 (in lattice units), exists such that Tc = for gH > gHc- This 
kind of behavior could be interpreted as the colored counterpart of the Meissner effect in 
ordinary superconductors, when strong enough magnetic fields destroy the superconductive 
BCS vacuum [37]. Then we shall refer to this remarkable result as the reversible color 
Meissner effect. Once again we would like to stress that this effect is not related to the 
color superconductivity in cold dense quark matter. Indeed, we believe that our reversible 
color Meissner effect is deeply rooted in the non-perturbative color confining nature of the 
vacuum and could be a window open towards unraveling the true nature of the confining 
vacuum. 

So far we reported our results for the critical temperature Tc in units of Aiatt and for the 
critical strength of the abelian chromomagnetic background field in lattice units. However 
it is well known that asymptotic scaling could be affected by scaling violation effects due 
to finite size of the lattice. On the other hand, such as effects are strongly reduced in the 
scaling of physical quantities. So that it is useful to analyze our data in physical units. In 
a pure gauge theory this can be done in terms of the string tension a computed at zero 
temperature in correspondence of the value of the gauge coupling (3 = (3c- We do not need 
to directly compute the string tension, for we may use the following parameterization of 
the SU(3) string tension given by Edwards et al. (see eq. (4.4) in ref. [52]) 

{ay/^){g) = fsu(?){9^) (l + 0.273Ia2(5)- 0.01545 a4(g) +0.01975 d6(5))/0.01364. 



"^^^"/5C/(3)(5^(/? = 6))' ^-g' 



(3.9) 



for 5.6 < /? < 6.5; fsu{3) is defined in eq. (|3.5|) . The critical temperature in physical units 
is given by 

Tr 1 , , 

(3.10) 



Moreover, using eq. ( 2.10 ), the field strength is 



/qH / Anriext 

y / ext ^g.^^^ 



y^(^ y L,a{^c) 

Our data for Td \fa versus \J gH j \fa on a 64'^ x 8 lattice are displayed in fig. ^ It is worth 
to note that, consistently with our previous analysis, lattice data can be reproduced by the 
linear fit 

= = a — (3.12) 
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with 

TciO) 



a 



0.643(15) a = -0.245(9) . (3.13) 



It is noticeable that our determination for Tc{0)/^/a is consistent with the determinations 
Tc/\/o' = 0.640(15) obtained in the Hterature without external field [51]. Using eq. ( 3.12| ) 
the critical field can now be expressed in units of the string tension 



^'^^ 2.63 ±0.15. (3.14) 



/a 

Assuming ^/a = 420 MeV, eq. (|1|) gives for the critical field 



= (1-104 ± 0.063)GeV (3.15) 

corresponding to gHc = 6.26(2) x 10^^ Gauss. Recently, it has been suggested that strong 
magnetic fields of order 10^^ Gauss are naturally associated with the QCD scale [53]. 
Moreover, it is believed that large magnetic fields might be generated during cosmological 
phase transitions. So that, we see that our findings could imply interesting effects during 
the cosmological deconfinement transition, which are worthwhile to investigate. 

3.2 SU(2) 

We also studied the SU(2) lattice gauge theory in a constant abelian chromomagnetic field. 
Even in this theory the deconfinement temperature turns out to depend on the strength of 
the applied chromomagnetic field, as already discussed in sect. ||. 

We evaluated the critical coupling P*{Lt,ncxt) on a 64^ x 8 lattice versus the strength 
of the external chromomagnetic field, introduced on the lattice by constraining the links 
according to eq. ( ^.11 ). As in previous section the critical coupling has been found by 



locating the peak of the derivative of the free energy density with respect to the gauge cou- 
pling p. Figure ^ shows our analysis for Tc in units of Aj^tt versus the critical temperature 
aTc together with a linear extrapolation to the continuum. As one can ascertain there is 
evidence for a dependence of the critical temperature on the applied field strength. As in 
the case of SU(3) the critical temperature can be expressed in terms of a physical scale 
by using a parameterization for the SU(2) string tension obtained by means of a fit to the 
string tension data collected in Table 10 of ref. [51]. We interpolate the string tension data 
by using Chebyshev polynomials of the first kind up to order 6 (see fig. P). 

In fig. ^ Tc/ \fa is plotted against ^ gH j \fo- As in the SU(3) case discussed in previous 
section, we can try to fit the data by means of a linear law. Remarkably we found that the 
linear fit eq. ( |3.12| ) works quite well and we get 

T fO) 

' 0.710(13) a = -0.126(5). (3.16) 



The value obtained for Tc(0)/-^/a is in good agreement with the value y/a = 0.694(18), 
without external field, obtained in the literature [51]. 
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Figure 5: SU(2) in (3+1) dimensions. Tc/Aiatt versus aTc for three different values of the external 
field strength (rioxt = 2,3,5) on 64^ x Lt lattices (Lt = 6,7,8). Full points in correspondence of 
aTc — are the extrapolations to the continuum. Full square is the critical temperature for SU(2) 
lattice gauge theory without external field taken from rcf. [50] . 



Now we can estimate the critical field in string tension units that turns out to be 

= 5.33 ± 0.33 . (3.17) 



Note that the critical field y/gWd ^fa is about a factor 2 greater than the SU(3) critical 
value in eq. ( |3.14| ) . This is at variance of the effective approach within dual superconductor 
picture in ref. [54], where one gets for the dual critical magnetic field gHc/a = 1 for SU(2), 
while gHc/a = 3/4 for SU(3). 

3.3 U(l) 



In sections p.l|and 3.2 we reported our results indicating a dependence of the deconfinement 



temperature on the strength of a constant abelian chromomagnetic background field. The 
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Figure 6: SU(2) string tension in (3+1) dimensions. Open circles are taken from Table 10 of 
ref. [51]. The solid line is our best fit with Chebyshev polynomials of the first kind up to order 6. 

main aim of this section is to find out if the effect we found is peculiar of non abelian gauge 
theories. To this purpose we consider four dimensional U(l) lattice gauge theory. 

It is known that, at zero temperature, U(l) lattice gauge theory undergoes a weak 
first order phase transition [55-57] from the confined phase to the Coulomb phase for 
(5 = 1.0111331(21) (using the standard Wilson action). We would like to seek a possible 
dependence of the confinement-Coulomb phase transition on the strength of an applied 
constant magnetic field. 

The quantity we have measured to locate the critical coupling is the derivative of the 
vacuum energy density (with respect to the gauge coupling) in presence of the background 
field (see sect. ^) 



where < U^y >»tcxt is the average plaquette evaluated with next 7^ and next = respec- 
tively. 

In fig. ^ we display the above quantity for three values of the constant abelian back- 
ground field, normalized to Sg^ti the derivative of the classical energy due to the external 
applied field 
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Figure 7: SU(2) in (3+1) dimensions. The critical temperature Tc estimated on a 64'^ x 8 lattice 



in units of the string tension, eq. (3.1C), versus the square root of the field strength \/gH in units 
of the string tension, eq. ( 3.11 ). Solid line is the linear fit eq. ( 3.12| ) On the zero vertical axis are 
represented the extrapolation of our data to zero value of the field (open circle) and the value for 
Tc{Q)/y/a (without external field) given in ref. [51] (full circle). The green full circle is the critical 



field in units of the string tension, eq. (3.17) 



The values of (3 corresponding to the peak in e'(/3, ^ext) for several values of the strength of 
the applied constant abelian field are displayed in fig. |9|. Our conclusion is that, contrary 
to non abelian lattice gauge theories, the critical coupling does not depend on the applied 
magnetic field strength. Analogous result was found in ref. [58] for (2+l)-dimensional 
compact QED. 



4. (2+1) dimensions 

Our numerical results for non abelian gauge theories SU(2) and SU(3) in (3+1) dimensions 
in presence of an abelian constant chromomagnetic background field lead us to conclude 
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that the deconfinement temperature depends on the strength of the apphed field, and 
eventually becomes zero for a critical value of the field strength. A natural question arises 
if this phenomenon, which is peculiar of non abelian gauge theories, continues to hold in 
(2+1) dimensions. To this purpose we consider here the non abelian SU(3) lattice gauge 
theory to be contrasted with the abelian U(l) lattice gauge theory at finite temperature. 

4.1 SU(3) 

In this section we focus on gauge systems in (2+1) dimensions. As is well known gauge 
theories in (2+1) dimensions possess a dimensionful coupling constant, namely has 
dimension of mass and so provides a physical scale. 

Non abelian gauge theories in (2+1) and (3+1) dimensions are sufficiently similar. 
Indeed, lattice simulations provide convincing evidence that (2+1) dimensional SU(N) 
gauge theories confine with a linear potential [59]. Moreover, at finite temperature there 
is a deconfinement transition [60]. 
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Figure 9: U(l) in (3+1) dimensions. The critical coupling (3c evaluated on a 64 x 16^ lattice 
versus the strength of the constant background magnetic field (in lattice units). The value at zero 
external field is the infinite volume extrapolation given in ref. [55]. The solid line represents the 
central value of /3c from ref. [55]. 



In (2+1) dimensions the chromomagnetic field H"^ is a (pseudo)scalar 

H'^ = ^e,,F,^^=F^,. (4.1) 
For SU(3) gauge theory a constant abelian chromomagnetic field can be obtained with 



exp(i2^^) O" 

exp(-iS^) 











1 



(4.2) 



As in the four dimensional case (see sect. p.3[ ) since we assume to have a lattice with toroidal 
geometry the field strength is quantized 

29H 27r_ 



a-— = —next , next integer . 
Z L\ 



(4.3) 
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Figure 10: SU(3) in (2+1) dimensions. The derivative of the free energy density with respect to 
the gauge coupUng f3, eq. ( ^.14| ), versus (3 for several values of the external field strength. Lattice 
size is L X 256 x Lt with two values of L = 256, 512 and temporal lattice size Lt = 4. Solid lines 



are the fits eq. (3.1). 



We computed the derivative of the free energy density eq. ( |2.5D on a L x 256 x 4 lattice, 
with L = 256,512 and several values of the external field strength parameterized by next- 



Our numerical results are reported in fig. |10|. We locate the critical coupling j3c as the 



position of the maximum of the derivative of the free energy density at given external field 
strength. As for SU(3) in (3+1) dimensions, the value of /3c depends on the field strength. 
Using the parameterization for the string tension given in eq. (C9) of ref. [59] 

^ ^ .X 4.1(1.7) 46.5(11.0) 

paV^ = 3.367(50) + ^ + ^ — ^ (4.4) 

we are able to estimate the critical temperature Tc in units of the string tension. We find 
that, as in (3+1) dimensions, Tc/y^ depends linearly on the applied field strength (see 
fig. 0). The hnear fit eq. ( p^ ) gives 

T (0) 

' - 1.073(87) Q = -0.193(76) , (4.5) 



(7 



that implies a critical field in string tension units \/ gild \fo = 5.5 it 3.7. Note that value 
for Tc(0)/y^ in the present work is in fair agreement with Td ^fo = 0.972(10) without 
external field obtained in ref. [60]. To check possible finite volume effects, we performed a 



lattice simulation with Lt = 8. The result, displayed in fig. 11, shows that within statistical 
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Figure 11: SU(3) in (2+1) dimensions. The critical temperature Tc estimated on 256^ x 4, 512^ x 4 
and 512 x 256 x 8 lattices in units of the string tension, eq. ( 3.10| ), versus the square root of the 
field strength \/ gH in units of the string tension, eq. (3.11). Open circles refer to Lt — 4, diamond 
to Lt — 8. Solid line is the linear fit eq. ( 3.12 ). In correspondence of \/gH / ^Ja ~ 0: full circle 
represents Tcj \fo at zero external field obtained by the linear extrapolation eq. ( 3.12 ), full square 
is the value given in ref. [59]. 



uncertainties our estimate of the critical temperature from the simulation with = 8 is 
in agreement with result at = 4. 

4.2 U(l) 

In a classical paper [61] Polyakov showed that compact quantum electrodynamics in (2+1) 
dimensions at zero temperature confines external charges for all values of the coupling. 
Moreover it is well ascertained that the confining mechanism is the condensation of mag- 
netic monopoles which gives rise to a linear confining potential and a non-zero string 
tension 

oa^ = exp[-7rV(0)/3] (4.6) 

VP 

where k is a constant, V{0) = 0.2527 [62] is the value of the lattice propagator at zero 
separation, and (3 = l/{ag^). 

At finite temperature it is well known that the gauge system undergoes a deconfinement 
transition which appears to be of the Kosterlitz-Thouless type [63]. We are interested in 
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Figure 12: U(l) in (2+1) dimensions. The derivative of the free energy density with respect to (3, 



eq. ( 3.18 ), versus f3 for several values of the strength of the constant magnetic background field on 
a 512 X 256 x 4 lattice. Solid lines are the fits to the data near each of the peaks using eq. (3.1). 
The inset is a magnification of the peak region. 



lattice U(l) gauge theory in an uniform external magnetic field 

C/f^*(x) = 1, 

\ (4.7) 
U2 (x) = cos{gHxi) + ism{gHxi) . 

We performed numerical simulations on 512 x 256 x 4 and 512 x 64 x 8 lattices. We measure 
the derivative of the free energy density with respect to the coupling /?. In fig. ^ we display 
the results for the 512 x 256 x 4 lattice. To determine the critical coupling /3c, we fitted 
the lattice data to eq. ( |3.lD . Contrary to the case of (2+1) and (3+1) non abelian lattice 
gauge theories, we do not find a dependence of the critical value of the coupling Pc on the 
magnetic field strength. Indeed we found that (temporal size Lf = 4) 

/3c(next = 5) = 1.694(17) 
/3c(next = 7) = 1.701(13) 
/3c(next = 9) = 1.716(11) 
/3c(next = 11) = 1.719(10)) 



By increasing the temporal size to Lj = 8 the critical coupling increases and is still inde- 
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Figure 13: U(l) in (2+1) dimensions. The derivative of the free energy density with respect to /3, 



eq. (3.18), versus (3 for several values of the strength of the constant magnetic background field on 



a 512 X 256 x 8 lattice. Solid lines are the fits to the data near each of the peaks using eq. (3.1) 



pendent of the external magnetic field strength (see fig. O). Indeed, we found 



/3c(next = 8) = 2.040(14) 

/3c(next = 9) = 2.051(15) (4.9) 
/3c(next = 10) = 2.054(12) 

Therefore we can conclude that even in (2+1) dimensional case the critical coupling does 
not depend on the strength of the external magnetic field as for U(l) lattice gauge theories 



in (3+1) dimensions (see sect. 3.3). 



5. Conclusions 



Let us conclude this paper by stressing our main results. We have investigated U(l), 
SU(2), and SU(3) pure gauge theories both in (3+1) and (2+1) dimensions in presence of 
an uniform (chromo)magnetic field. For non abelian gauge theories we found that there is 
a critical field gHc such that for gH > gHc the gauge systems are in the deconfined phase. 
Moreover, such an effect seems to be generic for non abelian gauge theories. On the other 
hand our numerical results for abelian gauge theories, where it is well established [32, 61] 
that confinement is due to monopole condensation, do not show any dependence of the 
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Figure 14: SU(3) in (3+1) dimensions. The derivative of the free energy density with respect to 
the gauge couphng /3, eq. ( 2.14 ), versus (3 at fixed external field strength (riext = 1) for a 64^ x 8 
lattice. Open symbols refer to T3 abelian chromomagnetic field, induced by lattice links in eq. ( |2.§| ). 
Full symbols refer to Ts abelian chromomagnetic field eq. (pj|). 



critical coupling from the strength of an external magnetic field. Therefore it seems very 
difficult to explain our reversible color Meissner effect in SU(2) and SU(3) gauge theories 
in terms of abelian color magnetic monopoles. Instead, the peculiar dependence of the de- 
confinement temperature on the strength of the abelian chromomagnetic field gH could be 
naturally explained if the vacuum behaved as an ordinary relativistic color superconductor, 
namely a condensate of color charged scalar fields whose mass is proportional to the inverse 
of the magnetic length. However, the chromomagnetic condensate cannot be uniform due 
to gauge invariance of the vacuum, which disorders the gauge system in such a way that 
there are not long range correlations. Consequently we can speculate that if the vacuum 
behaved as a non uniform chromomagnetic condensate, our reversible color Meissner effect 
could be easily explained, for strong enough chromomagnetic fields would force long range 
color correlations such that the gauge system gets deconfined. One might thus imagine the 
confining vacuum in non abelian gauge systems as a disordered chromomagnetic conden- 
sate which confines color charges due both to the presence of a mass gap and the absence 
of long range color correlations, as argued by R.P. Feynman in (2+1) dimensions [38]. 
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A. SU(3) Tg abelian chromomagnetic background field 



As is well known, in SU(3) there are two independent ways to realize a constant abelian 



chromomagnetic field. The first one, that we have considered in section |3.1| , is to take the 
abelian field directed along direction 3 in color space; the second one is to take the field along 
direction 8 in color space. In this Appendix we consider SU(3) lattice gauge theory in (3+1) 
dimensions in presence of a constant chromomagnetic background field along direction 8 
in SU(3) color space and along spatial direction 3. The continuum gauge potential and the 



corresponding lattice links are defined in eq. (2.7) and eq. ( |2.9| ) respectively. One should 
not expect a vastly different behavior in the two cases. Indeed we found that, even for 
a constant abelian chromomagnetic background field directed along color direction 8, the 
critical deconfinement temperature Tc depends on the strength of the applied field. Even 



more fig. [ij, where we display the derivative of the free energy density with respect to the 
gauge coupling f3 for ricxt = 1; 2 on a 64^ x 8 lattice, shows that the derivative of the free 
energy density in presence of the abelian chromomagnetic background field directed along 
color space direction 8 behaves like the case in which the field is directed along color space 
direction 3. Moreover the critical couplings at fixed field strength are consistent within 
statistical errors. 
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